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By using both the Lewis- Riesenfeld invariant theory and the invariant-related unitary transfor- 
mation formulation, the present paper obtains the exact solutions to the time-dependent supersym- 
metric two-level multiphoton Jaynes-Cummings model and the Chiao-Wu model that describes the 
propagation of a photon inside the optical fiber. On the basis of the fact that the two-level multipho- 
■ ton Jaynes-Cummings model possesses the supersymmetric structure, an invariant is constructed 

' in terms of the supersymmetric generators by working in the sub-Hilbert-space corresponding to 

a particular eigenvalue of the conserved supersymmetric generators ( i.e., the time-independent 
invariant). By constructing the effective Hamiltonian that describes the interaction of the pho- 
ton with the medium of the optical fiber, it is further verified that the particular solution to the 
Schrodinger equation is the eigenfunction of the second-quantized momentum operator of photons 
£^ ■ field. This, therefore, means that the explicit expression (rather than the hidden form that involves 

the chronological product ) for the time-evolution operator of wave function is obtained by means 
of the invariant theories. 
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I. INTRODUCTION 



The model that describes the interaction between a two-level atom and a quantized single-mode electromagnetic 
field is termed the Jaynes-Cummings (J-C) model [1], which can be applied to the study of many quantum effects such 
as the quantum collapses and revivals of the atomic inversion, photon antibunching, squeezing of the radiation field, 
inversionless light amplification, electromagnetic induced transparency [2-5], etc. By making use of the generalized 
invariant theory [6], we can obtain exact solutions and geometric phase factor of the two-level J-C model whose 
parameters of the Hamiltonian are totally time-dependent. Additionally, there exists another type of J-C model 
(so-called two-level multiphoton Jaynes-Cummings model) which possesses supersymmetric Lie-algebraic structure, 
fn this generalization of the J-C model, the atomic transitions are mediated by k photons [7,8]. Singh has shown that 
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■ this model can be used to study the multiple atom scattering of radiation and the multiphoton emission, absorption, 
and the laser processes [9]. Some authors introduced a supersymmetric unitary transformation to diagonalize the 
Hamiltonian of this J-C model and obtain the eigenfunctions of the stationary Schrodinger equation [10,11]. In 
the present paper, we generalize this method and obtain the exact solutions and the expression for the geometric 
phase factor of the totally time-dependent two-level multiphoton Jaynes-Cummings (TLMJ-C) model through the 
invariant-related unitary transformation formulation. 

The exact solutions and the geometric phase factor of the time-dependent spin model have been extensively inves- 
tigated by many authors. Bouchiat and Gibbons discussed the geometric phase for the spin-1 system [12]. Datta et 
aJ. found the exact solution for the spin- ^ system [13] by means of the classical Lewis- Riesenfeld theory and Mizrahi 
calculated A-A phase for the spin- ^ system [14] in a time-dependent magnetic field. The more systematic approach 
to obtaining the formally exact solution for the spin- j system was proposed by Gao et &1. [15] who made use of 
the Lewis- Riesenfeld quantum theory [16]. In this work, the generators of the Hamiltonian form the SU(2) algebra. 
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The dynamical algebraic structure of a type of general case of spin model such as the time-dependent L — S coupled 
system has been investigated and a set of SU(N) generators was constructed to linearize the Hamiltonian by Cen et 
al. [17]. Investigation in the direction of many-spin system has been completed by F. Yan et al. [18] who used a unitary 
transformation formulation which originated in the work of Gao in 1991 [6]. They discussed the time evolution of the 
Heisenberg spin system and obtained the formally exact solutions of the Schrodingcr equation in a time-dependent 
magnetic field. The relationship between the unitary transformations of the dynamics of quantum systems with 
time-dependent Hamiltonians and the gauge theories has been found by Montesinos et al. [19]. In particular, they 
showed that the nonrelativistic dynamics of spin- ^ particles in a magnetic field B % (t) can be formulated in a natural 
way as SU(2) gauge theory, with the magnetic field B l {t) playing the role of the gauge potential. This geometric 
interpretation provides a powerful method to find the exact solutions of the Schrodingcr equation. Although many 
investigations have been done in the direction of the time-dependent spin model, the subjects of all these works were 
the first-quantized cases. In this paper, we consider another time-dependent second-quantized spin model, namely, 
the generalized Chiao-Wu model [20], which describes the propagation of the photon inside the optical fiber. 

In 1984, Berry showed that the wavefunction in quantum adiabatic process would give rise to the geometric phase 
factor [21]. Afterwards, Chiao and Wu suggested their model, which is the first physical realization of Berry's 
geometric phase. It is well known that the geometric phase factor appears only in the systems with the time-dependent 
Hamiltonian or the Hamiltonian involving the evolution parameters. Both the classical geometric phase factor and 
the quantum geometric phase factor that is associated with vacuum, exist in the systems of the second- quantized 
field theory. In 1980'S, there was an argument concerning whether the geometric phase in the fiber experiment 
performed by Tomita and Chiao [22] belongs to the quantum or classical level [23-27]. Chiao- Wu'S theory [20] which 
is concerned with the polarization of the light propagating inside the noncoplanar optical fiber has no expressions for 
the Hamiltonian; moreover, their theory is first-quantization. From the point of view of us, it is consequently not 
applicable to investigating the geometric phase factor of quantum level. Only the system, of which whose Hamiltonian 
is second-quantized will present the geometric phase factor of quantum level. In what follows we show that the 
propagation of the photon in the fiber is essentially a second-quantized problem, i.e., the quantum-theory problem 
whose Hamiltonian is time-dependent. 

Note, however, that Berry's theory of the geometric phase proposed in 1984 is applicable only to the case of adiabatic 
approximation [21]. The invariant theory that is appropriate to treat the time-dependent systems was first proposed 
by Lewis and Riesenfeld (L-R) [16] in 1969. In 1991, Gao et al. generalized the L-R invariant theory and proposed 
the invariant-related unitary transformation formulation [6,28]. One of the advantages of this unitary transformation 
method is that it can transform the hidden form (which is related to the chronological product) of the time-evolution 
operator U(t) into the explicit expression. Many works have shown that the invariant-related unitary transformation 
approach is a powerful tool for investigating time-dependent systems and geometric phase factor [29-31]. 

II. THE INVARIANT THEORY AND THE INVARIANT-RELATED UNITARY TRANSFORMATION 

FORMULATION 

For the sake of illustrating the L-R invariant theory [16] easily, we consider a one-dimensional system whose 
Hamiltonian H(t) is time-dependent. According to the L-R invariant theory, a Hermitian operator I(t) is called 
invariant if it satisfies the following invariant equation ( in the unit % = 1 ) 



dl(t) 1 



+ 7 [I(t),H(t)}=0. 



(2.1) 



dt 



The eigenvalue equation of the time-dependent invariant \X n ,t) is given 



I{t)\X n ,t) 




(2.2) 



where 



= o. 



(2.3) 



The time-dependent Schrodingcr equation (in the unit % = 1) for the system is 



..0|*(t)) 



S 



= H{t) !*(*)) 



(2.4) 



dt 
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In terms of the L-R invariant theory, the particular solution |A„, t) of Eq. (2.4) differs from the eigcnfunction | X n , t) 
of the invariant I{t) only by a phase factor exp[i^ n (£)], then the general solution of the Schrodinger equation (2.4) 
can be written as 

l*(*)>a = EC n exp[Un(t)] \Xn,t) , (2.5) 
n 1 

where 

/■*/,, a ,\ , 

X n ,t ) dt , 



C n = (\ n ,t = 0\V(t = 0)) s . (2.6) 

(t)] \X n ,t) (n = 1, 2, • • •) are said to form a complete set of the solutions of Eq. (2.4). The statement 
outlined above is the basic content of the L-R invariant theory. 

A time-dependent unitary transformation operator can be constructed to transform J(i) into a time-independent 
invariant Iy = V^(t)I(t)V(t) [6,28] with 

Iv |A„) = A„ |A„) , (2.7) 

\X n ) =V^t)\X n ,t). (2.8) 



Under the unitary transformation V(t), the Hamiltonian H(t) is correspondingly changed into Hv{t) 

H v {t) = V\t)H{t)V{t)-V\t)i^. 



(2.9) 



In accordance with this unitary transformation method [6], it is very easy to verify that the particular solution |A„, t) sQ 
of the time-dependent Schrodinger equation associated with Hy (t) 

i dlXn dt t)s0 =H v (t)\X n ,t) s0 (2.10) 
is different from the eigenfunction |A„) of Iy only by the same phase factor exp[i<p n (t)] as that in Eq. (2.5), i.e., 

|A„,t) s0 = exp[^ n (i)]|A n ). (2.11) 
Substitution of |A n ,i) s0 of Eq. (2.10) into Eq. (2.11) yields 

|A n ) = H v (t) \X n ) , (2.12) 

which means that Hy(t) differs from Iy(t) only by a time-dependent multiplying c- number factor. It can be seen from 
Eq. (2.12) that the particular solution of Eq. (2.10) can be easily obtained by calculating the phase from Eq. (2.12). 
Thus, one is led to the conclusion that if the V(t), Iy, Hy(t) and the eigenfunction |A„) of Iy have been found, 
the problem of solving the complicated time-dependent Schrodinger equation (2.4) reduces to that of solving the 
much simplified equation (2.10). This paper obtains the exact solutions of the time-dependent Schrodinger equation 
describing TLMJ-C model and expression for its geometric phase factor by making use of this invariant-related unitary 
transformation method. 

III. THE EXACT SOLUTIONS OF THE TIME-DEPENDENT TLMJ-C MODEL 

The time-dependent Schrodinger equation and the Hamiltonian of the time-dependent TLMJ-C model under the 
rotating wave approximation are given by 
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H(t) = Lu(t)a^a + ^-a z + g{t){a^) k a_ + g*(t)a k a + (3.1) 

where and a are the creation and annihilation operators for the electromagnetic field, and obey the commutation 
relation [a, a*] = 1; a± and a z denote the two-level atom operators which satisfy the commutation relation [a z , a±] = 
±2(T-|- ; g(t) and g*(t) are the coupling coefficients and k is the photon number in each atom transition process; 
LVo(t) and ui{t) are respectively the transition frequency and the mode frequency. All the parameters in H{t) arc 
time-dependent. 

The supersymmetric structure can be found in the TLMJ-C model by defining the following supersymmetric trans- 
formation generators [10,11]: 

t fc-1 1 /ata+| \ > ( a k {a^) k \ 
* = °.+ 2 = { aat-|J' N ={ (at) * a * J ' 

Q = (at)^_=( (a ? )fe °), (/ (I;';;). (3.2) 

It is easily verified that (N, N , Q, Qt) form supersymmetric generators and have supersymmetric Lie algebra proper- 
ties, i.e., 

Q 2 = (Qt)2 =0j [Q\Q] = Ncr z , [n,n']=o, [N,Q} = Q, 
[JV,Q+]=-Qt, {Qt ;Q } = Ar ' ; {Q ;(Tz } = {Qt )(7z } =0, 

[Q,a z ]=2Q, [Q\a z ]=-2Q\ (Q* - Q) 2 = -N (3.3) 

where {} denotes the anticommuting bracket. By the aid of (3.2) and (3.3), the Hamiltonian (3.1) of the TLMJ-C 
model can be rewritten as 

H(t) = u(t)N + "W- S W az + g(t)Q + 5 *(t)Qt _ ^ (3.4) 

with S(t) = ku>(t) — u>o{t)- 

The equation which governs the time evolution of the TLMJ-C model is the time-dependent Schrodingcr equation 
(2.4). We will show the solvability of Eq. (2.4) by using the generalized invariant formulation in what follows. 

According to the invariant theory, we should first construct an invariant I(t). From the invariant equation (2.1) one 
can see that I{t) is the linear combination of N,<r z ,Q and Q' . However, it should be pointed out that the generalized 
invariant theory can only be applied to the study of the system for which there exists the quasi-algebra defined in 
Ref. [14]. It is easily seen from (3.3) that there is no such quasi-algebra for the TLMJ-C model. We generalize 
the method which has been used for finding the dynamical algebra 0(4) of the hydrogen atom to treat this type of 
time-dependent models. In the case of hydrogen, the dynamical algebra 0(4) was found by working in the sub-Hilbert- 
space corresponding to a particular eigenvalue of the Hamiltonian [32] . In this paper, we will show that a generalized 
quasi-algebra can also be found by working in a sub-Hilbert-space corresponding to a particular eigenvalue of N in 
the time-dependent TLMJ-C model . This generalized quasi-algebra enables one to obtain the complete set of exact 
solutions for the Schrodinger equation. It is easily verified that N commutes with H{t) and is a time- independent 
invariant according to Eq. (2.1). 

Use is made of a k {a)) k \m) = (m +, fc)! |m) and (a^) k a k \m + k)= (m +, fc)! \m + k), then one can arrive at 

W, '"L^vf, (3-5) 



to + k) J \\m + k) 

with A m = te±£H Thus we obtain the supersymmetric quasi-algebra (N,Q,Q\o z ) in the sub-Hilbert-space corre- 
sponding to the particular eigenvalue X m of N , by replacing the generator TV with A m in the commutation relations 
in (3.3), namely, 

[Q\Q] - \ m <T z , {Q*,Q} = \m, (Q t -0) 2 --A„ l . (3.6) 
In accordance with the invariant theory, the invariant I(t) is often of the form 

I(t) = c(t)g t + c* (t)Q + b(t)a z (3.7) 
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where c*(t) is the complex conjugation of c(t), and b(t) is real. Substitution of the expressions (3.4) and (3.7) for I{t) 
and H{t) into Eq. (2.1) leads to the following set of auxiliary equations 

c-^[cS + 2bg] = 0, c* + i [c*S + 2bg*} = 0, 

b+\\ m {c*g-cg*) = 0, (3.8) 

where dot denotes the time derivative. The three time-parameters c, c* and b in 7(f) are determined by these three 
auxiliary equations. 

Using the invariant-related unitary transformation method [6], we define the unitary transformation operator as 
follows 

V(t)=exp[a(t)Q-a*(t)Q^] (3.9) 

with a*(t) being the complex conjugation of a(t). With the help of the commutation relations (3.3), it can be found 
that, by the complicated and lengthy computations, if a(t) and a*(t) satisfy the following equations 

• /a *\ \i X m (ca*+c*a) „ i 

sin(4aa A m ) 2 = 1 t^, cos(4aa A m ) 2 = b, (3.10) 

(4aa*A m ) 2 

a time-independent invariant can be obtained as follows 

I v = VHt)I(t)V(t) = a z . (3.11) 

From Eq. (3.10), we substitute the time-dependent parameters A and 7 for c, c* and b in I{t) for simplicity and 
convenience, and the results are 

\ exp(i 7 ) t I cxpH 7 ) 
a = - r , a = ^ t . 

_ sinAcxp(i7) » _ sinAexp(-i7) 

C — J , C — J • ^O.i/J 

\ 2 \ 2 

Am Am 

Thus, the invariant J(i) in (3.7) can be rewritten 

I(t) = ^ r ^[cxp(«7)Q + cxp(-i7)Q t ] +cosAct z . (3.13) 
x 2 

Am 

In the meanwhile, under the unitary transformation (3.9), the Hamiltonian (3.4) can be transformed into 

H v (t) = V\t)H{t)V{t) - Vi(t)i^V{t) 

= ujN - ^ + {^(1 - cos A) + ^A m [gexp(-i7) + g* exp(z7)] sin A + 

+ ^-^cosA + |(1 -cosA)}cr z (3.14) 
where we use the Baker-Campbell-Hausdorff formula [33] 

V ^4t v ^ = wt L + h& L > L i + + h®m L > Q> L] > L] + ■ ■ ■ (3 - 15) 

with V(t) = cxp [L(t)]. The eigenstates of a z corresponding to the eigenvalue a = +1 and a = — 1 are ( Q ) and (J), and 
the eigenstate of N' is (|Jj™^) in terms of (3.5). From Eq. (2.6), (2.11), (2.12), we obtain two particular solutions of 
the time-dependent Schrodinger equation of the time-dependent TLMJ-C model which are written in the forms 

|* ro ,a=+i(*)> - cx V {-j\^ d , a=+1 (t) + ^ a=+1 (t)}dt}V(t) (3.16) 
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with 



and 



and 



with 



and 



<j>±*=+i(t ) = (m+ ^M*') + ^\i{g(t )exp[-i-y(t)]+g*(t) expl^t')}} sin \(t) 

-^icosA(i') (3.17) 



K*=+i(t) = 2^i[i_cosA(t')]; (3.18) 



|* m ,a=-i(*)> = cxp{i /"W=-i(*') + fll<7= _i(t')]dt'}V(t) f ° ) (3.19) 



0d,a=_i(t') = (m+ ^)w(t') + ^A|,{ff(t')exp[-i7(t')] + fl *(t')exp[*7(t')]}sinA(t') 

+ ^^cosA(i') (3.20) 



6 fl>tr =_i(t') = -^[1 - cosA(i')]. (3.21) 



These two particular solutions of the Schrodingcr equation (2.4) contain corresponding dynamical phase factor 
exp[i J 4>d,a{t )dt ] and the geometric phase factor exp[i J Q j> g , a (t )]dt with a = ±1. It should be noted that when 
the parameter A is taken to be a constant in the expression (3.18) and (3.21), the geometric phases in a cycle in the 
parameter space of the invariant I(t) can be rewritten as 

<j>a{T) — ^27r(l — cos A), (3.22) 

where 2n(l — cos A) is the solid angle over the parameter space of the invariant I(t) which represents the geometric 
meanings of the phase factor. Apparently, it can be seen that the former is dependent on the transition frequency u>o{t) 
and the mode frequency u>(t), and the coupling coefficients g(t) and g*(t) as well, whereas the latter is immediately 
independent of these frequency parameters and the coupling coefficients. 



IV. THE SECOND-QUANTIZED SPIN MODEL AND THE SOLUTION TO THE TIME-DEPENDENT 

SCHRODINGER EQUATION 

The exact solutions of the time-dependent spin model can be obtained by using the invariant-related unitary 
transformation method. The Hamiltonian of the spin model is generally of the form 

H(t)=c(t)-J, (4.1) 

where c(t) is the time-dependent arbitrary vector parameters and can be taken 

c(t) =co{t)[sm9(t)cos(p(t),sm6(t)sm<p(t),cos6(t)]. (4.2) 

With J± — Ji ± i J 2 satisfies the commuting relations [J3, J±] = ± J±, [J+, J_] = 2J 3 . Then the expression (4.1) for 
H(t) can be rewritten 

H(t) = co(t){^ sin 6 (t) exp[-i(p(t)]J+ + ^smO(t)exp[itp(t)]J- +cos9(t)J 3 }. (4.3) 

Note that the analytic solution to the time-dependent Schrodinger equation governing the spin model is easily 
obtained from the previous solutions to the time-dependent supersymmetric two-level multiphoton Jaynes-Cummings. 
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It follows from (3.1) and (3.3) that if when u>(t) and N are chosen to be and 1, respectively, and the supersymmetric 
operators, Q\ Q, and ^cr z , are respectively replaced with J + , J_ and J3, then the above supersymmetric Jaynes- 
Cummings is formally changed into the spin model, and it is therefore convenient to solve this time-dependent spin 
model. In other word, we can introduce an invariant-related unitary transformation operator V(t) as follows 

F(*) = exp[/?(i)J + -/r(t)J_] (4.4) 

where the time-dependent parameters are taken to be 

/?(*) = -^exp[-i 7 (t)], P*(t) = -Mexp[i 7 (t)]. (4.5) 

By making use of the Baker-Campbcll-Hausdorff formula [33] and V(t) in expression (4.4), or immediately from Eq. 
(3.14), one can obtain Hy(t) from H{t) 

Hv {t) = V\t)H{t)V{t)-V\t)i^ 

— { co [cos Acos# + sin Asin#cos( 7 — <p)] + 7 (1 — cos A)} J3. (4.6) 

Analogous to Eq. (3.11), the invariant 

I(t) = isinA(t)cxp[-z 7 (t)]J + + i sin X(t) exp[i 7 (i)] J_ + cosA(i)J 3 (4.7) 

is changed into the time-independent operator, J3, i.e., 

I v (t) = V\t)I(t)V(t) = J 3 , (4.8) 

under this unitary transformation V(t). From the two expressions (4.6) and (4.8), one can see that Hv(t) differs from 
Iv only by a time-dependent c-number factor. Thus with the help of (2.10), (2.11) and (2.12), it is easy to get the 
general solution of the time-dependent Schrodinger equation (2.4) 

l*(*)>a = £ °m eM-<t>m(t)]V(t) \m) (4.9) 

m * 

with the coefficients C m = (m, t = |*(0)) s . The phase <p m (t) = (pd,m(t) + 4> g , m {t) includes the dynamical phase 

4>d,m{t) = J*(m\VHt')H(t)V(t)\m)dt 
= m J Q co(t ){cos X(t ) cos 0(t ) + 

+ sinA(i')sin0(t')cos[ 7 (i') - ip(t')]}dt (4.10) 

and the geometric phase 

<f> g>m (t) = J (™\ - ^tf)^p- l m ) d *' = m /oT(0[l - cosA(t')]dt'. (4.11) 



In what follows, when the time-dependent arbitrary vector parameters of the Hamiltonian is chosen to be a special 
form, the time-dependent second-quantized spin model is reduced to the generalized Chiao-Wu model that describes 
the photon propagating inside the fiber. Consider a noncoplanarly curved optical fiber that is wound smoothly on 
a large enough diameter [22], a photon propagating inside the fiber is along it at each point at arbitrary time, one 

can draw a conclusion that the eigenvalue of the helicity • J of the photon is conserved in motion and its helicity 

operator ■ J is an invariant. In terms of the invariant equation (2.1), we construct an effective Hamiltonian H e f / (t) 
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which represents the interaction between the photon field and medium of the fiber. The effective Hamiltonian H e ff(i) 
is of the form 1 

k(t)x k (t) -, 

H eff( t ) = p J , ( 4 - 12 ) 

where J = Z + S is the total angular momentum operator of second-quantized electromagnetic field, which is given 

by 

L = - % -Y.S dk{[(k x Vg)at(fc, A)]a(fc, A) - [(fc x V s )a(k, A)]a+(fc, A)}, 

^ A 

S = - % -$ difef [at(fc, l)a(k, 2) - a){k, 2)a(k, 1) + a(fc, 2)a+(jfe, 1) - a(fc, l)at(fc, 2)]. (4.13) 
The infinitesimal rotation operator of motion of the photon in the fiber is given as follows 

U R = l-i$-J (4.14) 

with 

3= n ® X J®&, (4.15) 



where dot denotes the time rate of change of k(t). It follows from Eq. (4.14) that the photon state a, k(t)j ( a 
denotes the eigenvalue of helicity ) satisfies the following time-evolution equation 



d 
i — 



a, k(t) 



k(t)x k (t) f 



dt k 2 



■ J 



a,k(t)) (4.16) 



with H e ff(t) — fc (*)* 2 fc (*) . j being the effective Hamiltonian of photon inside the fiber. Apparently, together with 
the invariant I(t), the effective Hamiltonian H e ff(t) agrees with the invariant equation (2.1). From the form of the 
effective Hamiltonian H e ff(t) and the expression (4.13) for J, one can see that the problem of the rotation of the 

polarization plane is actually in analogy with that of the time-dependent second-quantized spin model, and fc ( t ) fc x a fc ( t ) 
may be considered a general magnetic field. 
Set the components of momentum of a photon 

^ = (sinA(t)cos7(t),sinA(f)sin7(f),cosA(t)) (4.17) 
k 



1 As far as we are concerned, the photon propagation problem can be ascribed to a time- dependent second-quantized spin 
model, where the effective (phenomenological) Hamiltonian is of a second-quantized form. Whereas in the previous researches 
[20-26], this problem was treated often by using classical Maxwell's Equations and first-quantized Schrodinger equation (and 
Berry's adiabatic geometric phase formula as well [21]). Although these investigations can be said to be somewhat outstandingly 
successful in both predicting and studying adiabatic geometric phases of photons in the fibre, here I still want to emphasize 
two points: for one thing, only by using the second-quantization formulation can we investigate the photon geometric phases 
at quantum level; for another, only when we consider the non-normal-product second-quantized Hamiltonian can it enable us 
to predict the existence of geometric phases at quantum-vacuum level. Tomita and Chiao may also agree to the above first 
point. They held the arguments [26] that although the geometric phases in the curved fibre can also be obtained by means 
of classical Maxwell's electrodynamics, they preferred to think of this phenomenon as originating at the quantum level, but 
surviving the correspondence-principle limit into the classical level. However, this point is not the main subject in the present 
paper, which will be further discussed elsewhere. Here, instead, we concentrate only on the second point, i.e., the geometric 
phases at quantum-vacuum level resulting from the zero-point radiation fields of vacuum, which has not been investigated in 
previous researches. 
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where the time-dependent parameters A(i) and j(t) denote the angle displacement of k(t) in the spherical polar 
coordinate system. In terms of the above process for getting the solution of the spin model, we use the invariant-related 
unitary transformation operator V(t) = exp[/3(i)J + — (3*(t)J-], which enables one to transform the time-dependent 
invariant I(t) into Iy(t) = Js- Under the unitary transformation V(t), one can obtain Hy(t) which is written as (4.6) 
from H{t) by making use of the Bakcr-Campbell-Hausdorff formula, where the time-dependent parameters 9 and ip 

represent the angle displacement of fc (*)* a fc M m ^ ne spherical polar coordinate, namely, 

k(t)x k (t) 

p = (sin 9 cos tp, sin 9 sin tp, cos 9). (4-18) 

By using Eq. (4.17), (4.18) and the invariant equation (2.1), two auxiliary equations can be derived 

7 sin 2 A = cos 9, A cos 7 — 7 cos A sin A sin 7 = sin 9 sin ip. (4.19) 
From Eqs. (4.19), it is shown that in Eq. (4.6) 

cos A cos 9 + sin A sin 9 cos(7 — ip) = 0, (4.20) 
thus, the expression (4.6) can be rewritten as 

ffv(*)=7(f)[l-cosA(t)]J 3 . (4.21) 
According to Eq. (4.11), the geometric phase of a photon whose eigenvalue of helicity is a can be expressed by 

M*) = {/o 7(*')[1 - cosA(t')]dt'} H J 3 \a) . (4.22) 

When we consider the adiabatic case, where A is time-independent time, t, the geometric phase <f) a in a cycle over the 
parameter space of the invariant I(t) may be written 

<t> a {T) = 2-7r(l — cos A), (4.23) 

where (<r| J3 \a) — 1 when the eigenvalue of helicity of the photon is taken to be a = 1. It should be noted that this 
result agrees with that derived from the Chiao-Wu theory [20]. 

By the aid of (4.9) and (4.22), it is easy to get the general solution of the time-dependent Schrodinger equation in 
the fiber experiment 

l*(*)> a - EC. expljMWit) \°) (4-24) 

a <■ 

with the coefficients C a = (a, t = |*(0)> s . 

It can be seen from the expression (4.20) and (4.21) that, in the noncoplanar optical fiber, the dynamical phase 
of photon due to the effective Hamiltonian vanishes, and its geometric phase is expressed by the expression (4.22). 
Since it is a second-quantized spin model, we do not take the normal product for the third component of the angular 
momentum J of the linear polarized photon. 



V. THE EXPLICIT EXPRESSION FOR THE TIME-EVOLUTION OPERATOR 

One of the applications of the invariant-related unitary transformation formulation is that it can change the hidden 
form of the time-evolution operator U(t) — P exp[i J Q H e ff(t )dt } of wave functions into an explicit expression, where 

P denotes the chronological product operator. Assume that the initial momentum k(t = 0) is parallel to the third 
component of the space coordinate, that is, 

fci =0, k 2 = 0, k 3 = k, (5.1) 

we have accordingly 

A(0) = 0, V(0) = 1. (5.2) 



9 



Use is made of the content of Sec. II which shows that 

U(t) = V(t)eMfJs), 

where 

/=|/o7(0[l-cosA(t')]dt'. 
At t = 0, the eigenstates of Jy corresponding to the eigenvalues cr = ±1 are of the forms 

e a^(fc) |0> , 



(5.3) 



(5.4) 



cr = +l,k(t = 0)) = —=[a){k,l)+ia){k,2)} |l 
' v2 

a = -1, fc(i = 0)) = ^[« f (fc, 1) - iat(fc, 2)] |0> = a{(k) |0) 



(5.5) 



where a R (k) and a^(fc) are the creation operators of the right- and left- rotation single-mode polarized photons, and 
k denotes the momentum of photon, and |0) represents the initial vacuum state at t = 0. Then cr, k(t)j at arbitrary 
time is given 



a, k(t)) = U(t) a, k(t = 0) ) = V(t) exp(/J 3 )4 |0> . 



In the following, we show that 



cr, fc(t)y is the eigenstate of the momentum operator 

p l i=J2f dkk^a' ! (k,X)a(k,\). 
x 



of photons field. Prior to this, some useful commuting relations are spread out as follows: 

[ip 2 + Pi, J+] = [ip2 -pi,J-] = 0, 

[1P2 - Pi, J+] = [iP2 + Pi, J-] = 2P3, 

[J3,Pn] = i(&lnP2 - 6 2 pPl), [P3, J±] = ip2±Pi, 
[Pm' j ±] = iS 3nP2 - ifauPs ± i(iSi»P3 - iS 3fl pi). 

Thus, the following formula can be calculated by using the above expressions 

p„tf(t)4(fc) |0) = V(t)\p^,eMfJ3)]aUk) |0) + V(i)exp(/J 3 )p„4(fc) |0) 
+ ]p ll ,V{t)]exp(fJ 3 )4{k)\0). 

First we compute exp(/ J 3 )]. By using the formulae (5.8) one will arrive at 



exp(-/J 3 )fvexp(/J 3 ) = i 



e-f — e f 



{S 2f iPi - 5niP2) + 



e f -|- e -/ _ 2 
+- 7; {SinPi + S 2 ^p 2 ) + P», 



then 



[p M ,exp(/J 3 )] = i 



exp(/J 3 )((5 2M pi - <5i M p 2 ) + 



+- 75 exp(/J 3 )((5i M pi + <5 2A1 p 2 ). 



Use is made of the expressions (5.5), we get 

Pn=J2f dkk ll a\k,\)a(k,\) = J dkk^{a R {k)a R (k) + a ] L (k)a L (k)} , 



(5.6) 



(5.7) 



(5.8) 



(5.9) 



(5.10) 



(5.11) 



then p^4(fc) |0) = fe Ai aJ.(fe) |0) . In terms of the initial conditions Eq.(5.2), the following eigenvalue equations of 
momentum can be obtained 
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p ia t(fc) |0) = fc sin A(0) cos 7(0)0* (fc) |0) = 0, 
p 2 al(k) |0) = fc sin A(0) sin 7(0)4 (fc) |0) = 0, 

p 3 aUk) |0) = fccos 7 (0)4(fc) |0> = kaUk) |0) . (5.12) 
By making use of the expressions (5.11) and (5.12), the first term on the right-hand side in (5.9) is evidently 

l/(*)[^,exp(/J 3 )]at (k) |0> = (5.13) 

and the second term is expressed by 

V(t)exp(fJ 3 )p^aUk) |0) - k^V(t) exp(/J 3 )o* (*0 |0) . (5.14) 
Next we calculate V(t)]. Using the commuting relations (5.8), the complicated calculation yields 
V^p^V = exp{-[/3(i)J+ - /3*(t)J_]}^exp[/3(t)J+ - /3*(t)J_] 

= + ^(A^ 3 + <^C) + ^(1 - cos A)(A M C - 4/?/3^3^3) (5.15) 

where 

C = /3*(pi-ijto) + /3(pi+ip2). (5.16) 

[p M , V(t)] can be derived from Eq. (5.15). Based on this, one can derive the third term on the right-hand side in (5.9), 
which may be rewritten as 

feV, V] cxp(/J 3 )4(fc) |0> = ^A^exp(/J 3 )p 3 4(*0 |0) 

-(1 - cosA)^Fexp(/J 3 )p 3 at (k) |0) . (5.17) 
According to Eq.(5.12) and (5.16), we can arrive at 

p„U{t)al(k) |0) = [k(S lf , cos 7 + <5 2aj sin 7 ) sin A - k{l - cos A)<5 3aj + k^U^a^k) |0) . (5.18) 

By combining Eq. (5.18) with (5.12), the eigenvalue equation of momentum of photons field can be rewritten in 
the form 

P ^U(t)aUk) |0) = k^(t)U(t)aUk) |0) (5.19) 

where fci(f) = fc sin A cos 7, k 2 (t) — fcsinAsin 7 , fc 3 (t) = fccosA. We thus show that U(t)a^(k) |0) is truly the eigenstate 
of momentum operator of the second-quantization photons field, and the eigenvalue is the time-dependent momentum 
fc M (t) of the single-mode photon. 

In accordance with the L-R invariant theory, the wave function a, k(t)\ = U(t)a^(k) |0) of the single-mode photon 



is also the eigenstate of the invariant I(t) = ^fi- ■ J. The significance of the invariant-related unitary transformation 
method will be illustrated in the following, by using the explicit expression for the time-evolution operator. The 
applications of the following commuting relations 

[J h V] = sinAcos 7 VJ 3 + COS7(1 ~ C ° SA V (/?*J-+/3J + ), 

A 

r -r ,„ . > . T . T sin 7 ( 1 cos A) , ^ ^ T , 

[J 2 V] = smAsm 7 W 3 + ^— '-V({3* J_ + /3J+), 

A 

[J 3 ,V] = (cos A - 1)W 3 + ^V(/TJ_ +/3J+) (5.20) 



enable one to get 



• J, V] exp(/J 3 ) = J 3 - V 1 ^- ■ Jexp(/J 3 ) 

= [/(i)J 3 -y/(i)exp(/J 3 ). (5.21) 
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Since 



M£ • JU(t)aUk) |0) - • J,exp(/J 3 )]at(fc) |0) + 



+ J,l/]exp(/J 3 )<4(fc)|0} + 

+U(t)^- Jat(fc)|0), (5.22) 
we apply the initial conditions (5.2) and (5.12) to Eq.(5.22) and obtain 

^ • JU(t)at (fc) |0) = aU(t)aUk) |0) (5.23) 

with the eigenvalue of • J being er = ±1 corresponding to the right- and left- rotation linear polarized photons, 
respectively. 

From what has been discussed above, one can draw a conclusion that one of the advantages of the invariant-related 
unitary transformation method is transforming the evolution operator U(t) of hidden form to the explicit expression 
(5.3). 

VI. CONCLUDING REMARKS 

We construct an invariant in the sub-Hilbert-space corresponding to a particular eigenvalue of the time- independent 
invariant N and get the exact solutions of the time-dependent TLMJ-C model by making use of the invariant-related 
unitary transformation formulation. In view of the above calculation, we can see that this unitary transformation 
formulation has some useful applications, for instance, it can solve the time-dependent systems and treat the geometric 
phase factor, and obtain the explicit expressions, instead of the hidden form, for the evolution operator of the wave 
functions. 

Since the three- level two-mode Jaynes-Cummings model plays an important role in Quantum Optics, the supersym- 
metric structure and the exact solutions of the time-dependent three-level two-mode multiphoton J-C model deserves 
further investigations by the formalism suggested in the present paper. 

We construct an effective Hamiltonian in this paper since the helicity of the photon field is a pseudo scalar and 
cannot be regarded as the Hamiltonian. We transform the problem of motion of the photon in the optical fiber into 
that of the time-dependent quantum spin model. Effective Hamiltonian is obtained by making use of the invariant 
equation (2.1), rather than through analyzing the electromagnetic interaction between the photons and the medium 
of the optical fiber. 

The invariant-related unitary transformation formulation is an effective method for treating the geometric phase 
factor [34,35]. This formulation replaces eigenstates of the time-dependent invariants by that of the time- independent 
invariants through the unitary transformation. It uses the invariant-related unitary transformation and obtain the 
explicit expression for the time-evolution operator, instead of the formal solution that is related to the chronological 
product. 



Note added: In addition, we also treat the non- normal-order spin operators and consider the potential effects (e.g., 
quantum-vacuum geometric phases) of quantum fluctuation fields arising in a curved optical fibre. The quantum- 
vacuum geometric phase, which is of physical interest, can be deducted by using the operator normal product, and 
the doubt of validity and universality for the normal-normal procedure applied to time-dependent quantum systems 
is thus proposed. Our brief history of investigating photon geometric phases in the fibre is as follows: in April 2000, 
Gao and I began to consider the non-cyclic non-adiabatic geometric phases of photon fields in the curved fibre 
based on a second-quantized spin model. In May 2000, Gao first proposed the concept of quantum-vacuum 
geometric phases. The existence problem of quantum- vacuum geometric phases is strongly relevant to whether 
the second quantization in spin model is adopted or not. Since it is in connection with properties of quantum 
electromagnetic vacuum and, moreover, this geometric phase is related close to the topological and global features of 
time evolution of vacuum-fluctuation fields, we think this concept is of essential significance and therefore deserves 
detailed investigations. From then on, these problems gained our attention and we tried to investigate this topological 
phases at quantum- vacuum level. 
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Since quantum-vacuum geometric phases has an important connection with vacuum energies, these experimental 
realizations may be relevant to the validity problem of normal-product procedure in the time-dependent quantum 
field theory (TDQFT), i.e., we also aim to re-examine the normal-product procedure in some extensions. If the 
quantum-vacuum geometric phases is proved present experimentally, then it is reasonably believed that it is not 
suitable for us to remove vacuum fluctuation energies and infinite charge density just by using the old formulation, 
e.g., re-defining the vacuum background energies and electric charges by utilizing the normal-product procedure, since 
in this re-definition, some potential physically interesting vacuum effects may also be removed theoretically. We think 
only for the time-independent quantum field systems can we use safely the normal-product procedure without any 
fear of introducing any new problems other than those which quantum field theory had encountered before. However, 
for the time- dependent quantum field systems, (e.g., photon fields propagating inside a helically curved fibre and 
quantum fields in an expanding universe), the physically interesting vacuum effects will unfortunately be deducted by 
the second-quantization normal-order formulation. So, the normal-order technique may therefore not be applicable to 
the time- dependent quantum fields. To the best of our knowledge, in the literature, this normal order problem in the 
time-dependent quantum field theory gets less attention and interests than it deserves. To test our above theoretical 
viewpoints, we hope the quantum-vacuum geometric phases of photons in the curved fibre would be investigated 
experimentally in the near future. 

Unfortunately, the left-handed polarized light due to vacuum fluctuation is often accompanied by the zero-point 
right-handed polarized light and their total quantum-vacuum geometric phases is therefore vanishing. So, it is not 
easy for physicists to investigate experimentally the quantum-vacuum geometric phases. This, therefore, means that 
our above theoretical remarks as to whether the normal-product procedure is valid or not for the time-dependent 
quantum field theory (TDQFT) cannot be examined experimentally. During the last three years, I tried my best but 
unfortunately failed to suggest an excellent idea of experimental realization of quantum-vacuum geometric phases. 
We conclude that it seems not quite satisfactory to test the quantum-vacuum geometric phases by using the optical 
fibre that is made of isotropic media, inhomogeneous media (e.g., photonic crystals 2 ), left-handed media (a kind of 
artificial composite metamaterial with negative refractive index), uniaxial (biaxial) crystals or chiral materials. Is 
it truly extremely difficult to realize such a goal? It is found finally that perhaps in the fibre composed of some 
anisotropic media (such as gyrotropic materials, including also gyroelectric and gyromagnetic media, where both elec- 
tric permittivity and magnetic permeability are respectively the tensors) the quantum- vacuum geometric phases may 
be achieved test experimentally. In these gyrotropic media, only one of the LRH polarized lights can be propagated 
without being absorbed by media. This result holds also for the zero-point electromagnetic fields. It is well known 
that people can manipulate vacuum so as to alter the zero-point mode structures of vacuum, which has been illus- 
trated in photonic crystals and Casimir's effect (additionally, the space between two parallel mirrors, cavity in cavity 
QED, etc.). If, for example, in some certain gyrotropic media one of the LRH polarized lights, say, the left-handed 
polarized light, dissipates due to the medium absorption and only the right-handed light is allowed to be propagated 
(in the meanwhile the mode structure of vacuum in these anisotropic media also alters correspondingly), then the 
quantum-vacuum geometric phase of right-handed polarized light can be easily tested in the fibre fabricated from 
these gyrotropic media. 



In order to illustrate our brief history (April, 2000 ~ March, 2003) of investigating geometric phases of photons in 
the noncoplanarly curved fibre, an outline is given as follows: 



2 Photonic crystals are artificial materials patterned with a periodicity in dielectric constant, which can create a range of 
forbidden frequencies called a photonic band gap. Such dielectric structure of crystals offers the possibility of molding the flow 
of light (including the zero-point electromagnetic fields of vacuum). It is believed that, in the similar fashion, this effect (i.e., 
modifying the mode structures of vacuum electromagnetic fields) may also take place in gyrotropic media. This point will be 
applied to the discussion which follows. 
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Chiao- Wu's model (photons moving inside a helically curved fibre) 



Berry's phase formula _ 



cyclic adiabatic geometric phase 4>f\T) — 2ira(l — cos A) Lcwis-Riesenfeld invariant theory _ 



non-cyclic non-adiabatic geometric phase c0\t) = crjj^jjt )[1 - cos Aft )]dt } second-quantized spin model ^ 



quantal geometric phases <j>^(i) — (tir — n L){jl j(t )[1 — cos X(t )]dt } 



under non-normal order 



quantum-vacuum phases <j)J=±i (t) = ±|{Jo j(t )[1 — cosA(i )]dt } 



validity problem of normal order in TDQFT 



experimental test is required I unfortunately, 0^ vacuum ) ft) + ^ vac ™) ft) = n by using gyrotropic-medium fibre 
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